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Abstract. This is an exposition of some of the aspects of quan- 
tum computation and quantum information that have connections 
with operator theory. After a brief introduction, we discuss quan- 
tum algorithms. We outline basic properties of quantum channels, 
or equivalently, completely positive trace preserving maps. The 
main theorems for quantum error detection and correction are pre- 
sented and we conclude with a description of a particular passive 
method of quantum error correction. 



1. Introduction 

The study of the underlying mathematics for quantum computation 
and quantum information is quickly becoming an interesting area of 
research |61j. While these fields promise far reaching applications [12, 
1411 160j , there are still many theoretical and experimental issues that 
must be overcome, and many involve deep mathematical problems. 
The main goal of this paper is to provide a primer on some of the 
basic aspects of quantum computing for researchers with interests in 
operator theory or operator algebras. However, we note that the only 
prerequisite for reading this article is a strong background in linear 
algebra. 

This work should not be regarded as an extensive introduction to 
the subject. Indeed, the reader with knowledge of quantum computing 
will surely have complaints about the selection of material presented. 
Moreover, we do not consider here the increasingly diverse fields of 
mathematics for which there are connections with quantum comput- 
ing (algebraic geometry, Fourier analysis, group theory, number theory, 
operator algebras, etc). Rather, our intention is to give a brief intro- 
duction and prove some specific results with the hope that this paper 
will help stimulate interest within the operator community. 
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The paper is organized as follows. The next section (§2) contains 
a discussion of the basic notions, notation and nomenclature used in 
quantum computing. In §3 we give a brief introduction to quantum 
algorithms by describing some elementary examples and presenting a 
simple algorithm (Deutsch's algorithm [19$ 20 ) that demonstrates the 
power of quantum computation. In §4 we outline the mathematical for- 
malism for the evolution of information inside quantum systems. This 
is provided by quantum channels, which are represented by completely 
positive trace preserving maps |15L I51L 1521 I62L I63j . The penulti- 
mate section (§5) includes a discussion of quantum error correction 
methods. We present the fundamental theorems for quantum error de- 
tection and correction in the 'standard model' |28l 1471 |64j. In the 
final section (§6) we describe a specific method of quantum error pre- 
vention |23L 1261 132L 14 7|, 1481 15 7|, I82j called the 'noiseless subsystem via 
noise commutant' method. Finally, we have included a large collection 
of references as an attempt to give the interested reader an entrance 
point into the quantum information literature. 



2. Quantum Computing Basics 

Let H be a (complex) Hilbert space. We shall use the Dirac nota- 
tion for vectors and vector duals in TC: A typical vector in 7i will be 
written as a 'kef and the linear functional on 7i determined by 
this vector is written as the 'bra' Notice that the products of a bra 
and a ket yield the inner product, (ipiWfa) , and a rank one operator, 
1^2) (^1 1- I n particular, given \ip) e TC, the rank one projection of H 
onto the subspace {X\ip) : A G C} is written \tp)(ip\. Further let B(7i) 
be the collection of operators which act on 7i. We will use the physics 
convention W for the adjoint of an operator U. 

The study of operators on Hilbert space is central to the theory of 
quantum mechanics. For instance, consider the following formulation 
of the postulates of quantum mechanics |61L 77J : 

(i) To every closed quantum system there is an associated Hilbert 
space 7i. The state of the system at any given time is described by a 
unit vector in 7i, or equivalently by a rank one projection 
When the state of the system is not completely known, it is repre- 
sented by a density operator p on 7i, which is a positive operator with 
trace equal to one. (This is the quantum analogue of a probability 
distribution.) 

{%%) The notion of evolution in a closed quantum system is described 
by unitary transformations. That is, there is a unitary operator U 
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on the system Hilbert space such that the corresponding evolution is 
described by the conjugation map p \— > UpW. 

(Hi) A measurement of a quantum system on TC is described by a set 
of operators M k , 1 < k < r, such that 

r 

E M l M * = 1 

The measurement is projective if each of the M k is a projection, and thus 
the M k have mutually orthogonal ranges. (A 'classical measurement' 
arises when all the projections are rank one.) The index k refers to the 
possible measurement outcomes in an experiment. If the state of the 
system is before the experiment, then the probability that event k 
occurs is given by p(k) = (if)\MlM k \ip) . Notice that {p(k)} determines 
a probability distribution. 

(iv) Given Hilbert spaces Tli, . . . ,H, m associated with m quantum 
systems, there is a composite quantum system on the Hilbert space 
Hi (g> • • -®H m . In particular, if the states of the individual systems are 
then the state of the composite system is given by \ipi)<S> - ■ •<S>\ip m )- 

The Hilbert spaces of interest in quantum information theory are of 
dimension N = d n for some positive integers n > 1 and d > 2. (It 
is generally thought that extensions to infinite dimensional space will 
be necessary in the future, but the current focus is mainly on finite 
dimensional aspects.) For brevity we shall focus on the d = 2 cases. 
Thus we let 7ijv be the Hilbert space of dimension 2 n given by the n-fold 
tensor product Hn = C 2 <8> . . . <8> C 2 = (C 2 )®". We will drop reference 
to N when convenient. Let (|0), |1)} be a fixed orthonormal basis for 
2-dimensional Hilbert space H 2 = C 2 . These vectors will correspond to 
the classical base states in a given two level quantum system; such as 
the ground and excited states of an electron in an atom, 'spin-up' and 
'spin-down' of an electron, two polarizations of a photon of light, etc. 
We shall make use of the abbreviated form from quantum mechanics 
for the associated standard orthonormal basis for 7i 2 ™ = (C 2 ) 0n ~ C 2 " . 
For instance, the basis for 7i 4 is given by {\ij) : i,j G Z 2 }, where \ij) 
is the vector tensor product \ij) = \i)\j) = \i) ® \ j). 

A quantum bit of information, or a 'qubit', is given by a unit vector 
= a|0) + in 7i 2 - The cases a = or b = correspond to the 
classical states, and otherwise \ip) is said to be in a superposition of the 
states |0) and 1 1) . A 'qudit' is a unit vector in C d . A vector state \ifi) in 
7Y/v is said to be entangled if it cannot be written as a tensor product of 
states from its component systems, so that \ip) does not decompose as 
IV') — I0i)l02) for some vectors i = 1,2. As an example, consider 
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\ip) = in 7^4, this is a vector from the so-called 'EPR pairs' or 

'Bell states'. Roughly speaking, the notion of decoherence'm a quantum 
system corresponds to the vanishing of off-diagonal entries in matrices 
associated with the system as it evolves. 

A number of specific unitary matrices arise in the discussions below. 
The Pauli matrices are given by 

x=r° iv y f ° - n ° 



i oj ' \i o y ' \o -i 

Let I2 be the 2x2 identity matrix. We regard these as matrix repre- 
sentations for operators acting on the given basis for Ti.2- In the n-qubit 
case (so N = 2 n ) we may consider all 'single qubit unitary gates' deter- 
mined by the Pauli matrices. This is the set of all unitaries {X^, Y^, : 
1 < k < n} where X x = X ® jf (n-1) , X 2 = % ® X ® lf n ~ 2 \ etc. Fur- 
ther let Ucn be the 'controlled-NOT gate' on TC4. This is the unitary 
which acts on 7i 4 by Ucn '■ ^ N) I (i + j) mod 2) . The CNOT gate 
has natural extensions {ujj^ '■ I < k ^ I < n} to unitary gates on Hn, 

where the fcth and Zth tensor slots act, respectively, as the control and 

(1 2) 

target qubits. For instance, with this notation Ucn = U C ' N . Note that 
Ucn only acts on the fcth and /th qubits in Ti^ = (C 2 )®". The set 

{Xk, Yk, Zk, Ucn '■ 1 < k ^ I < n} forms a set of universal quantum 
gates for T~Cn, meaning that this set generates the set oi N x N unitary 
matrices U(N) as a group (up to complex phases). We shall also make 
use of the Hadamard gate 

and the spin-| Pauli matrices = 1/2K for k — x,y, z. 



3. Quantum Algorithms 

Simply put, a quantum algorithm consists of an ensemble of initial 
states p which evolves under a unitary matrix U to a final density ma- 
trix UpW . The famous factoring algorithm of Shor [67, 70J and search 
algorithm of Grover |31| have received tremendous attention over the 
last decade. As a proper treatment of these algorithms is beyond the 
scope of this article, in this note we shall present the Deutsch algo- 
rithm [19] (and its generalization, the Deutsch- Josza algorithm [20J) 
since it is easily accessible and gives a good illustration of the power of 
quantum computation. In doing so, we shall give a description of two 
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fundamental classes of quantum algorithms: the simulation of a clas- 
sical function on a quantum computer and the algorithm for quantum 
parallel computation. 

Before continuing, let us illustrate a simple example of how an op- 
eration such as addition may be performed with a quantum algorithm. 
This will also allow us to establish some notation for the discussion 
which follows. We shall identify the standard basis vectors \i\ ■ ■ -i n ) 
for 7ijv with the integers {0, 1, . . . , 2 n — 1} via binary expansions. Then 
we may define a unitary U on Hn^'Hn by U \x)\y) = \x)\x®y) where 
the addition x © y is modulo N. The corresponding quantum algo- 
rithm implements addition modulo N. (Note that the CNOT gate is 
obtained in the case N = 2.) 

Towards the Deutsch algorithm, let H n = H® n be the n-fold tensor 
product of the Hadamard gate acting on Hn- Observe that H n \0)® n is 
the uniform superposition 

VZ x=0 

Fix positive integers k,m > 1. Let / : Z™ -> be a function. 
Consider the space 'H m ,,k = 7^2™ ® Ti.2 k - Then TC mj k has basis \x)\y) = 
\x)®\y) where x G Z™ and y G Z|. Define a unitary map U / : H m ,k — > 
^m,fc by 

C/j : \x)\y) i-> |x)|y©/(x)). 
For a given x G Z™, notice that the action of Uf is to permute the basis 
vectors : y G Z^}. 

Note 3.1. Observe that Uf\x)\0) = \x)\f(x)) for all x. Thus Uf simu- 
lates f on a quantum computer, and in this sense any classical function 
can be performed on a quantum computer. 

Next we compute 

2 m -l 

^((tf m io>n®ion = u f (-=j2^®w) 

2 m — 1 

= -i ® |0>) 

Hence an application of if m ® H 2 fc followed by [// yields a simultaneous 
parallel computation of / on all possible values of x. The corresponding 
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'circuit-gate' diagram for quantum parallelism for / is given below. In 
such a diagram, the 'circuits' correspond to states from the component 
systems (in this case |O) 0m and \0)® k from the systems 7i 2 m and H. 2 k 
respectively). The 'gates', drawn as boxes, indicate unitary operators 
applied as the system evolves with a left-to-right convention in the 
diagram (so here H m <g> l 2 k is applied first, then Uf is applied). 




|0> 



Let / : Z 2 — > Z 2 be a function. Then call / constant if /(0) = /(l) 
and balanced if /(0) 7^ /(l). The problem addressed by Deutsch's al- 
gorithm is the following: Given / : Z 2 — > Z 2 , determine if it is constant 
or balanced. Classically an answer to this problem requires two eval- 
uations of /. On the other hand, Deutsch's algorithm shows how to 
do this with one quantum operation. The circuit-gate diagram for the 
algorithm is given below. 



|0> 





H 




Uf 




H 












H 















The initial state is given by 1 01) = |0) ® |1). The first stage of the 
algorithm evolves this state to 



(H <g> H)(\0) ®\1)) = (H\0)) ® (H\l)) = |+) ® 



where |+) = and 
defined above yields 



-jtf-- Observe that the action of Uf as 



Uf(\x)®\-)) = (-l) f{x) \x)®\-) for xeZ 2 
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Thus, after the second stage of the algorithm the system has evolved 
to 

0>(l+>®|-» = ^0>((|o) + |i»®|-» 

^ (-i) /(0) |QK(-i) /(1) |i) j s| ) 

{±|+} ® |— ) if / is constant 
±|— ) <8> |— ) if / is balanced 

Finally, we apply a Hadamard gate H to the first qubit; that is, we 
apply H <g> I2 to the full system. Thus the system evolves to 

f (H ® 1 2 )(±|+) ® |->) = ±|0) ® |-) if / is constant 
\ (#® 1 2 )(±|-) ® |-)) = ±|1) ® |-) if / is balanced 

In particular, if we measure the first qubit we get 

J ±|0) if / is constant 
1 ±|1) if / is balanced 

Note that there is no uncertainty in the result: If we measure the first 
qubit and obtain ±|0) (respectively ±|1)), then we know / is constant 
(respectively balanced) with probability 1. 

Remark 3.2. The Deutsch-Josza generalization |2()| yields a more 
dramatic result. Let / : — > Z 2 be a function. Define / to be 
constant if fix) = f(y) for x, y G Z™, and balanced if |/ _1 (0)| = 
2 m_1 = |/ _1 (1)|. Suppose we know / is either constant or balanced 
and that we wish to determine which one / is. Classically, this requires 
2»n-i _|_ evaluations to know with certainty what type / is. The 
Deutsch-Josza algorithm shows how this may be accomplished with 
a single operation on a quantum computer. The algorithm acts on 
7^2 m ® ^2 with initial state |0)® m ® |1). The circuit-gate diagram may 
be obtained by adjusting the Deutsch diagram with |0)® m ® |1) in place 
of |0) <S> |1), m circuits at the top instead of one, H m in place of H in 
the top circuits prior to Uf, and H m in place of the final H. Thus, 

|O) 0m ® |1> ^ (H m ® h)U f (H m ® F)(|O) 0m ® |1)). 

Note 3.3. As a starting point for more recent work on quantum algo- 
rithms we mention |11L I59L I78| . Also see [7| for an extensive mathe- 
matical introduction to the study of quantum algorithms. 
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4. Quantum Channels 

While evolution in a closed quantum system is unitary (see postulate 
(ii)), experimentally evolution occurs in an 'open quantum system'. In 
such a system evolution is described mathematically by completely 
positive trace preserving maps [6lJ. The physical motivation for this 
is discussed below. 

A (linear) map £ : B{7i) — > B(7i) where 7i is a Hilbert space is 
completely positive if for all k > 1 the 'ampliation map' 

EM : M k {B{H)) -> M k {B(H)) 

given by £^ k \(pij)) = (£(pij)) is a positive map. (We could also write 
£(*) = %® £.) This is a rather strong condition; for instance, the 
transpose map is the standard example of a positive map which is not 
completely positive. The study of completely positive maps has been 
an active area of research in both the quantum physics and operator 
theory communities for at least thirty years. (In fact, it seems that 
many general results on completely positive maps have been obtained 
in the two fields without knowledge of the other.) See the texts of 
Kraus and Paulsen |62L I63| for good treatments of the subject 
from the two perspectives. 

Thus, in the general setting quantum information evolves through an 
open quantum system via completely positive trace preserving maps. 

Definition 4.1. A quantum channel £ : B(7i) — > B(7i) is a map which 
is completely positive and trace preserving. 

Trace preservation for a channel is equivalent to requiring the preser- 
vation of probabilities as states evolve through a quantum system. A 
channel is positive since density operators must evolve to density oper- 
ators, and it is completely positive because this property must be pre- 
served when the initial system is tensored with other systems (as part 
of a composite system). Physically, an open system can be regarded as 
lying inside a larger closed quantum system where all evolution occurs 
in a unitary manner. Thus, evolution in the open system can be re- 
garded as a 'compression' of the unitary evolution on the larger closed 
system. The mathematical formalism for this physical description is 
provided by Stinespring's dilation theorem [72J . 

The following fundamental result for completely positive maps was 
proved independently by Choi [15j and Kraus |52j. We present Choi's 
short operator proof. 
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Theorem 4.2. Let £ : B(TCn) —> B(Hn) be a completely positive map. 
Then there are operators E k G B(Hn), 1 < < iV 2 , sitc/i t/iat 

(1) S(p) = J2 E kpEl for all p e B(H N ). 

k=i 

Proof. Let eij = be the matrix units associated with the stan- 

dard basis for Hn- Let R = £^ N \(eij)). This matrix is positive by 
the iV-positivity of £. (In fact, Choi proved that the positivity of 
R characterizes complete positivity of £ .) Consider a decomposition 

R = SfcLi \ a k)( a k\, where \a k ) G C^ 2 are (appropriately normalized) 
eigenvectors for R. Let {Pj : 1 < i < N} be the family of rank 
iV projections on C^ 2 which have mutually orthogonal ranges and 
satisfy PiRPj = £[p%j)- Then |a&) = Y^i=i Pi\ a k)- Define operators 
E k : C N -> C w by £ fc |i) = P<|a fc >, so that 

Hence, 

£(ey) = £(|*>01) = iW 3 , = $> fc |*>0'l4. 

fc=i 

and equation (JI} holds by the linearity of £. ■ 

The decomposition is referred to as the operator-sum representa- 
tion of £ in quantum information theory. The operators E k are called 
the noise operators or errors of the channel. Note that trace preserva- 
tion of a channel is equivalent to its noise operators E k satisfying 

A; 

Remark 4.3. Choi's proof of this theorem provided the impetus for 
a recent application of Leung |56j to 'quantum process tomography'. 
This is a procedure by which an unknown quantum channel can be 
fully recovered from experimental data. To accomplish this for a given 
channel, (JI} shows that it is enough to recover the noise operators 
Ei. As Leung observes, Choi's proof shows the only experimental data 
required to fully reconstruct the E{ are the output states £{\i)\j)). This 
fact is the core of the tomography procedure outlined in |56j. 

The following result, which we state without proof, shows precisely 
how different sets of noise operators for the same channel are related. 
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Proposition 4.4. Let {Ex, . . . , E r } and {E[, . . . , E' r } be the noise op- 
erators for channels £ and £' respectively. Then 8 = £' if and only if 
there is an r x r scalar unitary matrix U = (uij) such that E = UE' 
where E t = [Ex ■ ■ ■ E r ] and (£")* = [E[ - -- E' r ]. In other words, 

r 

E i =^u ij E' j for l<i<r. 
i=i 

Let A = Alg{Ei, E}} be the algebra generated by the Ei and E\. 
This is the set of polynomials in the Ei and E\. An application of the 
Cayley-Hamilton theorem on minimal polynomials from linear algebra 
shows that all such polynomials may be written as polynomials with 
degree below some global bound. In quantum computing, A is called 
the interaction algebra for the channel. It is f-closed by definition, 
hence it is a finite dimensional C*-algebra [31 H T\ \7'6\ . Observe that, 
as a direct consequence of the previous result, A can be seen to be a relic 
of the channel; in other words, it is independent of the choice of noise 
operators which satisfy (JIJ for the channel. This is most succinctly 
seen in the case of unital channels, see § El for details. 

Note 4.5. We mention that an interesting and highly active area of 
current research in quantum information theory revolves around the 
study of quantum channel capacities. Specifically, there are a num- 
ber of deep mathematical problems which are concerned with comput- 
ing the capacity of a quantum channel to carry classical or quantum 
information. The following references give a starting point into the 
literature |51 12T1 1541 1551 15m 1551 HOI 1341 1551 1551 l6~ffl I7S| . 

Note that the quantum measurement process (postulate (Hi)) nat- 
urally determines a quantum channel. Let us consider more specific 
examples. (The final section includes further examples.) 

Examples 4.6. (i) Let < p < 1 and define operators on Ti.2 by 
Ex = (y/T=p)h, E 2 = (y/p)X, 

with respect to {|0), |1)}. The bit flip channel S(p) = ExpE\ + E 2 pE\ 
flips the state |0) to |1) and vice versa with probability p. For in- 
stance, observe that the projection |0)(0| evolves to the superposition 
5(|0)(0|) = (l-p)|0)(0|+p|l)(l|. 

(ii) Let Ex — y and define a channel by S(p) = E x pE x + a x pa x + 
o~ypa y + o z pa z . Observe that S(p) = 4f for every density operator 
p G M.2- This property, distinct density matrices evolving in a channel 
to the same density matrix, has recently been exploited as part of a 
scheme for quantum cryptography [5J. 
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(Hi) Let < r < 1 and define operators on 7i 2 by 

El = (o vT=t) ' E " = (o o) ■ 

These noise operators define an amplitude damping channel S(p) = 
E\pE\ + E 2 pE\- Such channels characterize energy dissipation within 
a quantum system. 

(iv) Let Ui, . . . , Ud be unitaries which act on a common Hilbert space 
and let ri, . . . , rd be positive scalars such that ^ r« = 1. Then we may 
define a channel by £(p) = J2i=i r iUipUj . These are the prototypical 
examples of 'unital' (S(l) = H) channels (see §EJ). In fact, every unital 
channel on Ai 2 may be written as a convex sum of unitaries in this 
way. This is not the case in higher dimensions however. 

(v) The class of entanglement breaking channels was introduced in 
|35j and studied in [39J. These are quantum channels which can be 
written in the form 

£{P) = l^fc)(^fc|(0fc|p|0fc)> 

k 

for some vectors \ipk) and \<pk)- With this representation, trace preser- 
vation is equivalent to Y^k \4>k)(4>k\ = 1, as Tr(p|0 fc )(0 fc |) = (<f) k \p\<j) k ) . 
Such channels derive their name from the fact that for d > 1, a density 
operator £^ d \T) is never entangled, even if T was initially entangled. 

5. Quantum Error Correction 

In this section we present the central aspects of the 'standard model' 
for quantum error detection and correction. For an extensive introduc- 
tion to the subject we point the reader to the articles [28, 14 7|, I64j . 
The general error correction problem in quantum computing is much 
more delicate when compared to error correction in classical computing. 
The possible errors that can occur include all possible unitary matrices, 
whereas in classical computing the only errors are bit flips. Nonethe- 
less, methods have been (and are being) developed for quantum error 
correction. 

5.1. Error Detection. Let Ti be the Hilbert space for a given quan- 
tum system. Then a quantum code C on 7i is a subspace of 7i. Let Pq be 
the projection of 7i onto C. Then Pq and P^ describe a measurement 
of the system which can be used to determine if a given state \ip) G Ti 
belongs to the code. The basic idea of an error-detection scheme in this 
setting is to first prepare an initial state in C, for brevity let us restrict 
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ourselves to unit vectors in C. The state \%p)(ip\ is then transmit- 
ted through the quantum channel 8 of interest, evolving to £(\ip)(ip\). 
Finally, the measurement Pq, P^~ is performed on this final state. This 
motivates the following definition. 

Definition 5.1. Let C be a quantum code on 7i and let E be an error 
(noise) operator associated with a given quantum channel on 7i. Then 
C detects the error E if the states which are accepted after E acts are 
unchanged, up to a scaling factor. In other words, there is a scalar \e 
such that 

P c E\tP) = A B |V) for all \ip) G C. 

Observe that the set of error operators E which are detectable for 
a fixed quantum code C form a subspace of operators, or a so-called 
operator space [62J. There are a number of equivalent conditions for 
detectable errors. 

Theorem 5.2. Let C be a quantum code and let E be an error operator 
associated with a given quantum channel. Then the following conditions 
are equivalent: 

(i) E is detectable by C, with scaling factor A#. 
(ii) P C EP C = X E Pc- 

(Hi) ^)x\E\th) = X E (ipi\\^2) for all |^> G C,i = 1,2. 
(iv) For every pair of vectors \ipi) and \ip2) in C which are orthogo- 
nal, the vectors E\tpx) and \ip2) are orthogonal. 

Proof. We shall prove the implication (iv) =>- (Hi). The other direc- 
tions are either trivial or easy to see. We may clearly assume that 
dimC > 2. Let {|'0i)> \ip2), • • •} be an orthonormal basis for C. We first 
claim that (iv) implies 

(2) {A\E\A) = for all 

Indeed, to see this, fix i,j and define 

|+) = |^> + |^-> and |-) = |^> - |^>- 

Then by (iv) we have 

= (+\E\-) = (^\E\^) - (^\E\^). 

Thus define A^ = (ipi\E\if)i) , and note that this is independent of i. 
Now let \ip) = axl^t) + a 2 |^ 2 ) + ... and \<j>) = /3i|^i) + W2) + ... be 
vectors in C. Then by (iv) and (J2J) we have 

(1>\e\<i>) = J2^(^\ E \^) = J2 a MHE\i>i) = A B (VII0>, 

hi i 

and this completes the proof. ■ 
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Let us describe a matrix perspective for detectable errors. This can 
be seen through a simple example. 

Example 5.3. Let C be the quantum code given by 

C = span{|000), |111)} C H 8 . 

The (unnormalized) error operators for the n-qubit depolarizing chan- 
nel |61l I44j are (E = {I2™, Zi, . . . , Z n }. In the 3-qubit case consider 
the error operator E = Z±. Observe that E\000) = 1 000) whereas 
= — I 111). Thus if E was detectable by C we would have, 

X E = (000|£|000) = 1 and X E = (111|£|111) = -1. 

From this contradiction it follows that none of the depolarizing er- 
rors Zk are detectable by the code C. In fact, from Theorem 15.21 the 
detectable errors for C may be realized in matrix form as the operator 
space I ( x } 2 * ) : A G Cj, with respect to an ordered orthonormal basis 
for H, s of the form {|000), 1 1 1 1) , . . .}, and where the * entries indicate 
that any choice is admissible. 

More generally, the conditions of Theorem 15.21 show that the de- 
tectable errors for a given quantum code C form the subspace of oper- 
ators 

' ' M * N | : A G C 

where the matrix form is given with respect to the spatial decomposi- 
tion H = P C H © PqH. 

5.2. Error Correction. Let <S = {Ei} be a set of errors that act as 
noise operators for a given quantum channel. If C is a quantum code, 
then the basic error-correction problem for C is to determine when 
there is a decoding procedure for C such that all the errors in (B are 
corrected. The simplest possible case occurs when every error Ei is 
the multiple of a unitary operator and the subspaces E£ are mutually 
orthogonal. The obvious decoding procedure in this situation is to first 
make a projective measurement to determine which of the subspaces 
EiC a given state £ C has evolved to, then apply the inverse of the 
error operator E^ This particular case motivates the following general 
definition in the standard model for quantum error correction. 

Definition 5.4. Let £ be a quantum channel and let C be a quantum 
code with projection Pq. Then C is correctable for £ if there is a 
quantum channel TZ such that 

1lo£(p)=p 

for all p supported on C; that is, all p with p = P c pP c . 
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There are a number of useful characterizations of correctable codes. 
Note that in the proof of (Hi) =>- (i) below, the error correction op- 
eration 1Z is explicitly constructed and hence this gives a constructive 
approach for decoding within this error correction model. 

Theorem 5.5. Let £ be a quantum channel with errors <£ = {Ei} 
and let C be a quantum code with projection Pq. Then the following 
conditions are equivalent: 
(i) C is correctable for £ . 

(ii) The operators in the set = {EIE2 : Ei G (£} are detectable 
byC. 

(Hi) There are scalars A = such that 

(3) P c E\E P c = \ lJ P c for all 

(iv) There is a linear transformation Ei i— > E[ on (E such that the 
new error operators E[ satisfy the following properties: 

(a) The subspaces E'fi are mutually orthogonal, 

(b) The restriction of every E[ to C is proportional to a restric- 
tion to C of a unitary operator. 

Proof. Observe that any matrix A which satisfies must be positive 
since this equation may be written as a matrix product A^A = (XijPc) 
where A = [E\Pc -E/2-Pc • ■ • ] is a row matrix. Conditions (ii) and (Hi) 
are equivalent by definition. We shall prove (i) is necessary and suf- 
ficient for (Hi) and leave the connection with (iv) for the interested 
reader. 

For (i) =^> (iii), let C be a quantum code with code projection P c . 
Suppose £ is a quantum channel with errors {Ei} and that C is cor- 
rectable for £ via the error-correction operation TZ with noise operators 
{Rj}. Define a compressed channel by £c(p) = £(PcpPc)- Then by hy- 
pothesis TZ(£c(p)) = lZ(£(PcpPc)) = PcpPc- In particular, 

Y,RiEiP c pPcE\R) = PcpPc for all p. 

Thus by Proposition 14.41 there are scalars such that 
RkEiPc = a ki P c for all i, k. 

Hence, 

(4) P c E\R\R k E 3 P c = a ki a kj P c for all i,j,k. 

But 7Z preserves traces, so that R\,R k = 1, and thus when we sum 
(0) over k we find 

P c E]E 3 P c = XijPc for all 
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where A„ = ^2 k a ki a kj . 

Conversely, let us assume that 8 is a channel with errors {-Ej} and C is 
a code with projection Pq such that holds for a positive scalar matrix 
A = {Ajj}. Let U be a unitary operator such that D = WAU = (d k i) 
is diagonal (so that D = di&g(d kk )). Note that ^2 k d kk = 1 by trace 
preservation of E. By Proposition 14.41 the operators F k = ^ u^Ei also 
implement E. A simple computation shows that 

PcFl^Pc = d kl P c for all fc,Z. 
The polar decomposition of F k Pc gives 



F k P c = U k \IP c FlF k P c = ^d~ k U k P c 



k' 

for some unitary U k . Define projections P k = U k PcU\. Then 

yon ntn U[P C F^ F k P c Ul 

P l P k = P l l P k = JL— — ^ = for fc^Z, 

V 

and hence the ranges of the P k are mutually orthogonal. 

Without loss of generality assume that ^2 k P k = H (otherwise just 
add the projection onto the orthogonal complement and define U k = I). 
The candidate error-correction operation is defined by 

n{ P ) = Y,u ] k P k pP k u k . 

k 

Observe that for all p with p = PcpPc, 

V "fcfc 



— 8 k i v d kk Pcp — S k iV d kk p. 

Thus we have 

ft(£(p)) = J2 U l P k F iP F l p kU k = J2 6 kid kk p = P 

k jit P& 

for all p = PcpPc, and we have proved that (Hi) =>■(«). ■ 
Let us discuss a pair of illustrative applications of Theorem 15.51 

Examples 5.6. (i) With C = span{|000), | 111) } inside H 8 as in Ex- 
ample 15.31 it is easy to see that C satisfies (jSJ for the errors £ = 
{1, Xi, X 2 , X 3 }. Let P be the projection onto C, let Pi be the projec- 
tion onto the subspace X\C = span{|100), 1 1 1 ) } and similarly define 
projections P 2 and P 3 . Then the correction operation given by the 
proof above is K = {P , X 1 P 1 , X 2 P 2 , X 3 P 3 }. 
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(ii) Shor's 9-qubit code |69| is defined by two orthonormal vectors 
in 7i 2 9 given by 

(|000) + |111))(|000) + |111))(|000) + |111>) 



\0l) 



2V2 

(|000) - |111))(|000) - |111))(|000) - 1 1 n>) 



2V2 

Let C = span{|0i), |1l)}- Given 1 < k < 9, the code C is correctable 
for the errors {X^, Z{\ as (jHJ) is satisfied for this triple. Let P x ^ be 
the projection onto the subspace X&C and similarly define projections 
Py^k and P Zj fe. Then the correction operation given by the theorem is 

71 = {XkP x ,k, YkPy,k, ZkP z ,k}- 

Since the Pauli matrices, together with the identity operator H2, form 
a linear basis for M. 2 that is closed under multiplication up to scalar 
multiples, it follows that C is correctable for any set of errors (E which 
act on one of the nine possible qubits. (In fact the Shor code also 
corrects for errors on multiple qubits |61j.) 

Note 5.7. The quantum error correction conditions of Theorems 15.21 
and !5.5l were established independently by Bennett, DiVincenzo, Smolin 
and Wooters [9| and Knill and Laflamme [50J . As a collection of entry 
point references for particular methods of quantum error correction we 
mention [1311^150113511351153115511711 155] . 

6. Noiseless Subsystems via the Noise Commutant 

In this section we describe a specific method of passive quantum 
error correction, by which we mean no active intervention is required 
after information is encoded. The basic idea in a 'noiseless subsystem 
method' of error correction, classical or quantum, is to encode infor- 
mation on subsystems of the system of interest in such a way that it 
remains immune to the effects of the channel that the information is 
evolving through. Let £ : B(H) — > B(7i) be a quantum channel with 
noise operators {Ei} and interaction algebra A = AlgjE'j, Ej}. Recall 
that £ is unital if £ (1) = £\ E t E\ = 1. The noise commutant for £ is 
the f- algebra 

A' = {p G B(H) :pA = Ap for A G A} 

= {peB(H):[p,E t } = 0=[p,Ej} for i = l,...,n}. 

The procedure in the noiseless subsystem via noise commutant method 
of quantum error correction f^l^Bll^lliTlliglfSTllgTl^ is to 

use the structure of A' to produce noiseless subsystems (which are also 
called 'decoherence-free subspaces' in special cases, see Remark I6.3J) . 
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The illustrations of this method that appear in the literature all 
involve unital channels. The following discussion shows why this is the 
case. Let Fix(£) = {p G B(TC) : £{p) = p}. Observe that Fix(£) is 
a f-closed subspace of B(7i). Now consider a unital channel £. Let p 
belong to A'. Then pEi = Eip for all i and hence 

£{p) = Y,EipE\ = p£{t)=p. 

i 

Thus p belongs to Fix(£). The converse inclusion was proved indepen- 
dently in [1 3 and [55] (see also [32] for another proof): The set Fix(£) 
coincides with the noise commutant in the case of a unital channel. In 
fact, it is not hard to show that the identity Fix(£) = A' characterizes 
unital channels. Thus, by building on the proofs from I32L 155] we 
may state the following. 

Theorem 6.1. Let £ be a quantum channel. Then Fix(£) = A' if and 
only if £ is unital. Moreover, in this case A is equal to the algebra Aq 
generated by the Ei. In other words, the operators E\ belong to Aq. 

Note 6.2. Notice how the algebra A = (A)" = Fix(£)' is a relic of the 
channel in the unital case (see the discussion after Proposition 14.41) . 

Remark 6.3. In the case of unital channels, Theorem 16. II shows that 
the noiseless subsystem via noise commutant method can be presented 
as follows: Let £ be a unital quantum channel. Then the interaction 
algebra A is generated by the Ei and A' = Fix(£) is a finite dimensional 
C*-algebra. As such, A' is unitarily equivalent to a unique direct sum of 
ampliated full matrix algebras, A' = @k{^m k ®M-n k )- Hence, a density 
operator p that encodes the initial states of a quantum system will be 
immune to the noise of the channel as it evolves through, £(p) = p, 
provided that p is initially prepared on one of the ampliated matrix 
blocks l mfc (g> M. nk inside the noise commutant A' = Fix(£). The 
decoherence-free subspace method may be regarded as the special case 
of this method that occurs when matrix blocks l mfe ®M. nk with m& = 1 
are utilized. 

In the case of a general quantum channel, however, the full struc- 
ture of A 1 cannot be used for error correction. This can be seen most 
dramatically in the following simple case. 

Proposition 6.4. Let £ be a completely positive map such that As = 
£(1) is not invertible. Let Pg be the projection onto the subspace Tig = 
Ran(/Lj). Suppose that Ei = PgEiPs for all i. Then is non-zero 
and every operator p G Biji.^) belongs to A' and satisfies £(p) = 0. 
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Proof. The non-invertibility of Ag implies TCg is non-zero. Let p > 
belong to B(Hg), by which we mean p acts on 7i with p = P^pP^. 
Then p trivially belongs to A' since pEi = = Eip for all i. Further, 
£(p) = £(P±p) = 0. ' ■ 



Remark 6.5. It would be interesting to know if the noise commutant 
can be used to produce noiseless subsystems for classes of non-unital 
channels. For instance, a natural generalization of unital channels is 
the class of channels for which the identity evolves to a multiple of a 
projection. Notice that if £(ln) — m P f° r some projection P, then m 
divides the dimension of 7i by trace preservation of £. A simple exam- 
ple of this phenomena is given by the channel £ with noise operators 
Ai = \0)(i\ for 1 < % < dim7Y = d. Trace preservation of £ may be 
readily verified, and in this case 

d d 

8(1*) = J>4 = £(|0><i|)(|i><0|) = rf |0)(0|- 

8=1 1=1 

Note 6.6. There are a number of other quantum error correction 
schemes that have been investigated and some of these will be of inter- 
est to operator researchers (See [11 1251125113511531157115011511 82]). 

We finish by considering noiseless subsystems for some special cases 
of unital channels. 

Examples 6.7. (i) Let < p < 1 and let Ei,E 2 be operators on 7i 2 
defined on the standard basis by Ei = — p)l 2 and E 2 = (y/p)Z. 
Then Ei and E 2 are the noise operators for a unital channel £ on 
M.2 which is a variant on the bit flip channel discussed earlier. The 
quantum operation corresponding to this channel is equivalent to the 
phase flip or phase damping operation on single qubits |61j. It is so 
named because, for instance, = (1 — +p\— )( — | and 

hence £ flips the phase of |+)(+| and |— )( — | with probability p. It is 
easy to see in this case that 

Fix($)=^' = {E 1 ,E 2 }' 

^ :a,6ecj ^C1©C1, 

and hence this channel has no non-trivial noiseless subsystems. 

(ii) Let < p < 1 and let Ei,E 2 be operators on TC4 defined on 
the standard basis by E\ = ^1 — p(l 2 ® 1 2 ) and E 2 = ^JpiZ ® Z). 
These noise operators determine a unital channel £ on M.^ which can 
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be regarded as an ampliation of the phase flip channel. Compute 

Fix($)=A' = {E U E 2 }' 

[fan a 14 \ 1 



Thus A is unitarily equivalent to the direct sum A = M.2 © M-2 and 
there is a pair of 2-dimensional noiseless subsystems. 

(Hi) More generally, let E\ = Ui,...,Ed = Ud be unitaries on a 
Hilbert space 7i. If we are given scalars Aj such that Yli M 2 = 1' 
then the operators \{Ei are the noise operators for a unital channel 
£ with Fix(£ ) = A' = {Ui, . . . , Ud}'- In this case the automatic self- 
adjointness of A can be seen directly through an application of the 
Cayley-Hamilton theorem. (The algebra generated by any unitary U 
can be seen to include W via minimal polynomial considerations.) 

(iv) An important special case of such channels is the class of 'col- 
lective rotation channels' p ITOllSmiSffl 1521 1351 13511751 ITffl I5HII5T] . 
The n-qubit example has noise operators given by weighted expo- 
nentiations of the operators Jk = Y^n=i^k f° r ^ = x iVi z i where 
Jfc = <Jk® (l2)^ n ~ 1 \ etc. There is an abundance of noiseless sub- 
systems for these channels and they can be computed directly |35J or 
through combinatorial techniques discussed below. 

(v ) A generalization of the collective rotation class is presented in 
|42j and noiseless subsystems for this class of 'universal collective rota- 
tion (ucr) channels' are computed using Young tableaux combinatorics. 
For each pair of positive integers d, n > 2 there is a family of such chan- 
nels with the case d = 2 yielding the class from (iv). It is proved in |42j 
that every ucr-channel possesses an abundance of noiseless subsystems. 
In fact, it is shown that the noise commutant for every channel in this 
class (for fixed d and n) contains the algebra A^ = Alg{7r(cr) : a G S n } 
where S n is the symmetric group on n letters and tt : S n — * U(d n ) is 
the unitary representation of S n on Tid^ given by 



for all hi, . . . , h n £ Hd and a £ S n . In particular, the Young tableaux 
machine can be used to explicitly compute the structure of A^ and 
identify noiseless subsystems for the corresponding ucr-channel. Fur- 
ther, a recent paper of Bacon, et al [4 J , has shown that the change of 
basis transformation from the standard basis to the Young tableaux 
basis can be implemented efficiently with a quantum algorithm. 
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